Solitonic approach to the dimerization problem in correlated one-dimensional systems 
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Using exact diagonalizations we consider self-consistently the lattice distortions in odd Peierls- 
Hubbard and spin-Peierls periodic rings in the adiabatic harmonic approximation. From the tails 
of the inherent spin soliton the dimerization dao of regular even rings is found by extrapolations to 
infinite ring lengths. Considering a wide region of electron-electron onsite interaction values U > 
compared with the band width 4to at intermediately strong electron-phonon interaction g, known 
relationships obtained by other methods are reproduced and/or refined within one unified approach: 
such as the maximum of doo at U~3to for t/~0.5 and its shift to zero for g— ><? o f»0.7. The hyperbolic 
tangent shape of the spin soliton is retained for any U and g ~0.6. In the spin-Peierls limit the doo 
are found to be in agreement with results of DMRG computations. 
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There is a longstanding debate on the interplay of 
the electron-electron (el- el) interaction and the electron- 
phonon (el-ph) interaction in conducting polymers like 
trans-(CH) x . Among a large amount of papers we refer 
here only to Refs. 1-6 devoted to the discussion of the 
origin of the observed dimerization d in the framework 
of the (extended) Peierls- Hubbard model (PHM). The 
PHM is regarded as the minimal microscopic model for 
conducting polymers. Special features of the interplay 
of on-site correlation U and off-diagonal el-ph interac- 
tion have been pointed out first in Ref. 1. Employing 
the Gutzwiller approximation (GA), it was shown that 
for weak and intermediate el-ph interaction strength, the 
dimerization d passes through a maximum near U ~4io 
when U is increased after which it is suddenly suppressed. 
The geminal approachO (GEA) shows a smooth decrease 
of d with increasing U in the opposite strongly corre- 
lated limit. According. tn_extrapolations based on exact 
diagonalizations (ED)ErEm the enhancement of d due to 
U predicted by the GA near the maximum is overesti- 
mated. In Refs. 2,3 the infinite chain limit c?oo calculated 
using the PHM was extrapolated from above (below) by 
2n-membered open chains (4ro+2- membered rings). In 
the highly correlated limit of the 1/2- filled band case the 
PHM can be mapped onto the antiferromagnetic spin- 
1/2 Heisenberg model (AFHM) and at low temperatures 
a spin-Peierls phase is expected. Such dimerized phases 
have been observed for CuGeOs (having Cu02 chains) 
and a'-NaV205. The CUO3 chains in Sr(Ca)2Cu03 are 
at the threshold to the AFHM-limitEI. In many cases of 
practical interest the actual d is very small and no self- 
consistent d ^0 is found in feasible short even rings. 

In the present paper we show that studies of odd rings 
of comparable lengths yield reliable estimates of lattice 



distortions in infinite rings at arbitrary strength of el-el 
correlations and reasonable strength of the el-ph inter- 
action. We exploit the generic property of neutral odd 
periodic rings that their ground state is given by a spin 
soliton. The lattice distortions for such a ring are shown 
schematically in Fig. 1. The bond in front of the soliton 
center, i.e. between sites N and 1, is a long bond for 
which aNi-aQ—2uo =d> holds ( a^i (ao) is the bond- 
length in the distorted (equidistant) state). For large 
N the region far from the soliton center tends to the 
regularly dimerized state. Hence, varying the model pa- 
rameters, insight into the behavior of d^ may be gained 
already at finite N from the study of d. For short rings 
(N = 3, 5, .), d(N) exceeds significantly d^. This is, the 
consequence of a strong first-order Jahn- Teller effectM In 
Refs. 8,9 the total energy E tot of odd AFHM rings has 
been studied for fixed geometries adopting rigid sharp 
soliton shapes. In our method all bonds are optimized 
to yield the minimum of the total energy E to t ■ The in- 
herent soliton exhibits a smooth shape. We illustrate 
our method considering the PHM and the AFHM, where 
comparison with other reliable approaches is possible. In 
particular, the density matrix renormalization group is 
used to check the extrapolated doo . 

For the electronic part H e \ of the total Hamiltonian 
H = H c \ + H\ at we adopt the one-band extended PHM 



+ Hi {Urii^n^i + V it i + inin i+ i) , 
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at half filling, where c\ s creates an electron with spin 
s=±l/2 at site i, ^i, s =cj s Cj S is the number operator, 
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and rii=J2s n is- We linearize the bond-length dependent 
transfer integral t and the intersite el-el interaction V 

k,i+l = - (to - IfVi) , Vi t i+i = V - r\Vi, (2) 

where Vi=u i+ i~Ui, itj is the displacement of the i^ 1 site 
relative to the undistorted state. In the adiabatic and 
harmonic approximation the lattice part H\ a t reads as 



(3) 



where K is the spring constant. Via the Hellmann- 
Feynman theorem we obtain N self-consistent Eqs. 



Kv t = A/N - 2 7 P M+1 + r)D iii+1 , 



(4) 



where A=J^ (27Pj j j + i — r)Di t i+i) expresses the fixed 
length constraint J^. = 0, -Pj.i+i being the bond or- 
der, and Di^ + i denotes the density-density correlators 
in the ground state | G) 
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(G | mni + i I G). 
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The strength of the el-ph interaction can be measured 
by the parameter g introduced asH 

g = 1 /^KT <\. (7) 
Note that sometimesB the related quantities defined as 
Ad-ph = 2g 2 /7r, or A S sh = 2A c i_ ph . (8) 



are used. Values of g ~ 0.4 to 0.5 typical for conducting 
polymers are regarded as weak to intermediate coupling 
constants. Following Refs. 5,6 we shall use hereafter the 
dimcnsionlcss dimcrization defined as 



d = 2u y / K/t < 1. 



(9) 



The parameter S frequently used to describe the modu- 
lation of the transfer integral i^+i = t (l + (— l)' l+1 S) in 
the dimerized statcEhlj is related to d by 5 = gd. In the 
limit U/to~S> 1 the low-energy physics of the 1/2-fiiled 
PHM ring (Eq. (Q)) can be described by the AFHM 



Ji,i+iSiS i+ i, with Jj^+i 
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(10) 



where the exchange integral is given to 2 nd order pertur- 
bation theory in to/U. Using Eq. (2) we find 

= Jo - j sp Vi = J (l + (-l) i+1 S sp ), (11) 

where r y sp =2j Jq / to and S sp =2gd characterizes the regu- 
lar spin-Peierls state for 77=0. We shall use the AFHM 



also out of the limit U^to for the case V=0 adopting 
an effective exchange integral J=J(U/t ) given by the 
relation Jij=(2/Tr)v sp tij and the spin velocity v sp taken 
from the Bethe-Apaatz solution for the equidistant infi- 
nite Hubbard ringed (hereafter ao=l, h=l) 



J, 
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where I n , (n = 0, 1) are modified Bessel functions. 

Using the Lanczos-method, the Hamiltonians (Eqs. 
([j], |Ic|)) have been diagonalized exactly for finite rings 
with periodic boundary conditions starting with a given 
set {v n ^}. The AFHM has been treated using the 
spinless fermion technique^ resulting in analogous self- 
consistent equations as for the PHM. Then the corre- 
sponding "ground-state" eigenvector | G) has been used 
to calculate the next set of lattice order parameters 
{v n iW\ using Eq. (|]). The iteration was continued un- 
til the maximal deviations of E tot and all u n ,(j) between 
two iteration steps j and j + 1 became smaller than the 
required accuracy of 10~ 7 to 10 -8 (see Refs. 3,12). Here 
rings composed of up to 7V=13 sites (PHM) and N=23 
(AFHM) have been studied. The ED computer limita- 
tions to rings where finite size effects are still important 
can be circumvented at least for even membered regu- 
larly dimerized AFHM rings and reasonable el-ph (sp-ph) 
interaction strength applying the density matrix renor- 
malization group techniqucO (DMRG) with typical dis- 
carded errors of the order 10~ 6 . 

At first we consider the one-particle Su-Schrieffer- 
Heeger model (SSH) (PHM: U,V=0) where very long odd 
rings can be treated numerically. The infinite even ring 
problem is reduced to a transcendental equation for doc: 
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k 2 ' 



with k = (.gdoo) 2 , (13) 



where K and E denote complete elliptical integrals. The 
results for up to TV = 601 sites are shown in Fig. 2. Start- 
ing from small N, d passes at first by a minimum at N 
~ 2£, the width of a soliton in an infinite ring. For weak 
el-ph coupling f > 1, cUn/doo « 0.86 at 2£/N w 0.92. 
Below that minimum all curves, for which N > 2£ holds, 
approach a nearly universal curve, pass a very sruall max- 
imum near 6£, and tend finally to doo from abovelij. Since 
for correlated problems only relatively short rings can be 
treated exactly, a strong nonmonotonous behavior could 
cause problems in extrapolating to doo . Fortunately, our 
calculations indicate that at least in the AFHM-limit the 
depth of the first minimum at finite N is strongly sup- 
pressed. 

Let us now consider how the on-site interaction U af- 
fects d. As shown in Fig. 3 starting from U = 0, d in- 
creases with U and has a maximum at t/ max ~ 3io, af- 
ter which d starts to decrease smoothly. Both behaviors 
are similar to the GA and GEA predictions, respectively. 
Quantitatively, however, we obtain J7 max (3)~3.12to for 
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g ^0.5. For 0.5< g <0.6, U max starts to decrease. Fi- 
nally, when g — > g c w0.7, U max — > 0. For compari- 
son we note that the GA's results are f/ max ~ 4/j for 
g < 5 C =0.76. The ED-results of Ref. 5 yield 5c =0.75± 
0.04 slightly above our result. Above g c , there is no 
enhancement of d due to U . Turning to larger rings, 
we first adopt an 1/N extrapolation (dashed curve) and 
arrive at a rough lower bound being crudest for small 
U. An improved bound is achieved connecting the ex- 
act U=Q point with the 1/N extrapolation upshifted to 
the AFHM limit (see below). To avoid an artificial min- 
imum, we omit the U=to point. Thus for g—0.5 we 
get 0.174<doo imax ^djv = i3 «0.2. Instead the GaI gives 
rfmax~0.31. For U > to, g <0.5 the dpHM — ^afhm from 
above. To compare pur d with the continuum model re- 
sult of Inagaki et aZJl3, we rewrite their dimerization as 
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where c=(4/7r) N /273«1.04 for k = 0.5 (see Eq. (10) and 
Refs. 16 -18). Applying Eq. ( [Ti] ) to intermediately cor- 
related cases, we adopt the effective exchange integral J 
defined by Eq. ( |i~2| ) and arrive at an analytic expression 
(the dashed-dotted curve in Fig. 3)E3. Surprisingly it ex- 
hibits similar shape and magnitude as the weak coupling 
(.9 < 0.5) ED-curves. In particular its [/ max /io=3.21 is 
close to 3.12 mentioned above. This suggests that even 
in the case of conducting polymers, being clearly outside 
the usual AFHM-regime, the dimerization is mainly gov- 
erned by the (always present) spin degrees of freedom 
and to less extent by the charge degrees of freedom. In 
the usual case U >2V, the V and its derivative rj (see Eq. 
(||)) enhance d. For U ^>to,V and ?/=0 one can replace 
U-> U-V. 

For long rings N > 2£, G?Ar-values close to doo can be 
expected. Then d = d(N) might be approximated by 



d(N) = d 



y? Al_ 

i=i 



exp 



—N 



(15) 



Note that in contrast with the general PHM case, d(N) 
for An and 4?i.+2 AFHM-rings can be described by one 
smooth curveEl. The d vcn — > ^oo from below, just oppo- 
site to odd rings. To be specific, we consider one typical 
example. We estimate for the upper curve shown in Fig. 
4, dodd.oo = 0.0765, Imax = 1, j4i,odd~0.5 and £ dd=5.26. 
The even ring curve tends from below to a slightly 
larger value d C vcn,oo=0.078 and ^4i, GVOn =-11.8. The expo- 
nent £ evon =1.85 differs significantly from £ dd- From the 
soliton shape (— l) n v n txd N tanh(ra/£jv), we deduced a± 
N=23, £n~5, djv«0.08, whereas the continuum modelO 
yields C=37r/j /(16J g 2 )=4.91. The fit of the even curve 
can be somewhat improved adopting l max — 4. Then 
with A\ = A 3 = and A 2 = -3.5, A± = -2040 one ar- 
rives at the same £ = 5.26 as in the odd case for l max = 1. 



Taking the DMRG-values for 7V=60, we conclude that 
the accuracy of the solitoniq-estimate of doo is ~ 2 to 
3 %. The continuum theoryl3, (Eq. (14)) predicts, for 
k=0.5, c?oo=0. 07203, a value slightly below our discrete 
results. According to our numerical finding we would 
recommend to use k «0.279. Fitting alternatively the 
curvature of d dd at large N by a parabola, one arrives 
at an extrapolated very shallow minimum at finite ring 
length (iV m i n w29 sites for the present parameters). Then 
the slightly smaller doo.odd compared with doo, even might 
be viewed as a hint for a tiny minimum at finite N generic 
for d oc id being the deepest in the SHH case (see Fig. 2). 
Raising g, the soliton becomes narrower. Thus at large 
g any minimum should be accessible by the ED. With 
increasing g the d a dd(l/N) curves become flatter. Small 
minima were detected for g=0.9, 0.85 at A r =7V m ; n =17, 
21, respectively. Anyhow, the l/TV-extrapolation of c? oc id 
(d e ven) from accessible N<N m in yields a lower (upper) 
bound of ■ 

To summarize, a novel approach to the dimerization 
problem of correlated ID-models has been presented. It 
is based on exact diagonalizations of odd ring Hamilto- 
nians combined with a self-consistent treatment of the 
classical lattice degrees of freedom. Known dependences 
of the bond alternation on the el-el and el-ph coupling 
strengths obtained by other approximations valid in dif- 
ferent parameter regimes have been reproduced and re- 
fined within one unified method. The 1/iV-extrapolation 
to the infinite rings gives a new lower bound for any cor- 
relation strength. With the aid of the Bethe-Ansatz so- 
lution for the spin velocity, even in the intermediate cou- 
pling regime a sizeable part of the dimerization can be 
described by an effective spin-Hamiltonian gaining thus 
new insights in the dimerization mechanism of conduct- 
ing polymers. The DMRG is found out as a valuable 
supplementary tool to our solitonic method. 
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FIG. 1 Schematical view of the lattice distortions in 1/2- 
filled odd rings. (Un)distorted sites are denoted by °(«). 



FIG. 2 Reduced dimerization in odd SSH rings vs. recipro- 
cal ring length 1/iV for various el-ph interactions g. The ring 
length N is given in units of the soliton width 2£=2 / '(gd<x>) ■ 



FIG. 3 Dimerization in the Peierls-Hubbard model vs. on- 
site energy U for the el-ph coupling g=0.5. In deriving the 
improved lower bound, Eqs. (13-15) have been used. 



FIG. 4 Size dependence of the dimerization d for even (o) 
and odd (•) periodic spin-Peierls rings. Even rings are treated 
by ED until A=22. The d for A=28 to 60 are obtained by 
the DMRG. The parameter set used 7 sp =0.4, J=l/3, K=l 
corresponds to 17=13, V=to=l,r]=0, and g=0.6 for the PHM. 
The full and the dashed curves are the fits by Eq. (Fl). 
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